Effect of magnetic impurity correlations on Josephson tunneling 
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The ordering trend of magnetic impurities at low temperature results in the frustration of the 
pair-breaking effect and induces a "recovery" of superconducting properties. We show that this 
effect manifests itself in the deviation of the Josephson current amplitude from the values obtained 
within the Ambegaokar-Baratoff and the Abrikosov-Gor'kov models. We consider both weak and 
strong-coupling cases. The theory is applied to describe the experimental data obtained for the 
low-T c superconductor SmRli4B4. We further predict a "recovery" effect of the Josephson current 
in high-temperature superconductors. 



I. INTRODUCTION 

The interaction between mobile charge-carriers and 
magnetic impurities is accompanied by a spin-flip pro- 
cess. In superconductors, this process leads to the 
pair-breaking effectErQ, that is, to the destruction of 
Cooper pairs. As a result, superconducting properties are 
strongly affected by the presence of magnetic impurities. 
For example, an increase of magnetic impurity concentra- 
tion riM leads to thej-decrease of the critical temperature 
T c and to gaplesnessaH, to a decrease of thexritical field 
-ff C 2U or the Josephson current amplitude JjM (at given 
temperature) . 

Most studies including the effect of magnetic impuri- 
ties consider the magnetic moments as independent from 
one another. Because of the conservation of the total 
spin (disregarding spin non-conserving couplings such as 
the dipole-dipolc interaction) scattering processes involve 
the spin-flip of both the charge carrier and the local- 
ized magnetic moment. However, in the region near 
T = one should take into account the correlations be- 
tween impurity spins. The ordering trend of the mag- 
netic moments frustrates the spin-flip process and de- 
creases, therefore, the pair breaking. This decrease leads 
to the "recovery" effect of superconductivity As shown 
by two of the authors in Ref. | the account of this effect 
explains the unusual temperature dependence of H c2 m 
high-temperature superconductors and leads to excellent 
agreement with all experimental data available so far. 

The interplay of superconductivity and ordering of 
impurity spins has also been studied in the context 
of ternary compounds containing rare fffth elements ( 
RERh 4 B 4 with RE =Sm, Ho, Er or hS u). In some of 
these materials (RE =Sm, Ho, Er) the impurity spin- 
spin correlations lead to a magnetic phase transition 
into an (anti)ferromagnetic state below a critical tem- 
perature T m . SmRh 4 B 4 is of special interest because 
it shows coexistence of superconductivity and antiferro- 
magnetism and displays an unusual temperature depen- 
dence of iJ C 2(T , )ErE3. Anomalies in the Josephson critical 
current have also been observed in this material E The 
experiment shows that the Josephson amplitude is en- 



hanced in the antiferromagnetic region as compared to 
the value extrapolated from the paramagnetic phase. 

In the following, we calculate the critical-current am- 
plitude J c of a Josephson junction in the presence of 
magnetic impurities and study the influence of the low- 
temperature ordering trend. Our analysis shows that in- 
cluding impurity spin-spin correlations leads to a sizeable 
and experimentally observable increase of J c at low tem- 
peratures with respect to the saturated value obtained 
when neglecting these correlations. It would be interest- 
ing to carry out experiments verifying our predictions. 
The experiment done on SmRIi4B4 is also discussed in 
the light of our theory. 



II. MAIN EQUATIONS 

A. The Josephson current 

Based on the method of thermodynamic Green's func- 
tions, one can obtain a set of equations for the supercon- 
ducting order parameter and the renormalization func- 
tion in the presence of magnetic impurities: 
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where D nn i = il 2 / [ft 2 + (uj n — o-i n ') 2 ] is the phonon 
Green's function, Q is a characteristic phonon frequency 
and A is the coupling constant. Here and in the fol- 
lowing we define A n = A(kj?,o; n ) and Z n = Z(k-p,u) n ) 
as the order parameter and renormalization function of 
the isotropic superconductor [u„ = (2n + 1)ttT with 
n = . . . , —1, 0, 1, . . .; k.F is the Fermi wave- vector] . The 
function u n that appears in Eqs. ([!]) and (|J) is solution 
of the equation (see Refs. |l]-§|): 
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T s ~ um is the spin-flip scattering amplitude (tim is 
the magnetic impurity concentration). In the absence of 
magnetic impurities (T s = 0) the previous expression re- 
duces to u n = u n /A n . Solving Eqs. (|l])-(||) for a given 
concentration of impurities (that is, for a given value of 
the parameter T s ) and for a given coupling constant A 
(or T c , see below) one obtains u n , which allows us to 
calculate the Josephson current. 

Consider a symmetrical Josephson junction and let us 
focus on the stationary Josephson effect. The amplitude 
J c of the Josephson critical current is given by (see, e.g., 
Ref. |§ 

n n 

where R is the normal state resistance of the barrier. 
This relation applies for both weak and strong-coupling 
superconductors. The current depends on the temper- 
ature and on the concentration of magnetic impurities 
through u n . 

In general, the sum over n cannot be evaluated ana- 
lytically. This is especially true for the strong-coupling 



case for which the dependence of u„, Eq. (^), on n ap- 
pears not only through u> n , but also through A„ and is 
thus non-trivial. However, there are two limiting cases, 
that are of interest for the present work, in which an an- 
alytical expression can be derived. Both are within the 
weak-coupling BCS theory (A„ = A). In the absence of 
magnetic impurities (EL = 0) one obtains the well-known 
Ambegaokar-BaratoffEJ result: 

where is Boltzmann's constant. In the numerical part 
we compare this temperature dependence with the one 
obtained for a superconductor containing (un) correlated 
magnetic impurities. Eq. (||) has two important features. 
The current saturates as T — > 0, and is linear near T c . 
These features remain when uncorrelated magnetic im- 
purities are added to the system. In the second limiting 
case we allow for the presence of magnetic impurities, but 
set T = 0. Then, 
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with f = r s /A. This expression is usefull for the de- 
termination of J c at T = in the weak-coupling case. For 
strong-coupling superconductors one has to use Eqs. (|l|) 
and (Q). The numerical solution near T = can then be 
obtained using the Poisson Formula. 



B. Spin Ordering of Magnetic Impurities 



The equations of the previous section describe the 
pair-breaking effect induced by the time-reversal non- 
invariant perturbation due to magnetic impurities. The 
theory accounts for such a perturbation through the spin- 
flip scattering amplitude T s = l/r s (see Eq. (||); t s is the 
relaxation time). in the usual picture of independent 
magnetic momentsEI, the amplitude T s is a temperature 
independent parameter. As already stressed in the intro- 
duction, the assumption of independent impurity spins 
is a good approximation at high temperatures and low 
impurity concentrations, but is not valid in the other 
cases. One example where the independent impurity- 
spin assumption is not valid is given by systems under- 
going a magnetic phase transition at a critical temper- 
ature T m . The observed deviations from the expected 
behaviour of superconducting properties near and below 
the (anti)ferromagnetic transition temperature of ternary 
compounds mentioned in the introduction point towards 
a change in the value of T s due to spin-spin correla- 
tions. We stress, however, that spin-spin correlations 



are important even if one does not observe a magnetic 
phase transition. For example, the formation of a spin 
"glass" is related to these correlations. The importance 
of magnetic impurity correlations increase with descreas- 
ing temperature. The anomal behaviour of the critical 
field H C 2 observed in the cuprates can be explained by 
such correlational Another manifestation of such &, "re- 
covery" effect is the decrease in microwave lossesE3 ob- 
served under application of an external magnetic field. 
In this case, the frustration of the spin-flip process is due 
to the external field. 

Because of the spin-spin correlations (the ordering 
trend of magnetic impurities) the spin-flip scattering am- 
plitude depends on temperature. This dependence T S (T) 
was studied by two of the authors in Ref. ^ and was ap- 
plied to the evaluation of H C 2- According to this study, 
r s (T) can be represented in the form: 
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where t = (T - T x )/6, T x = T m - ST with ST = a9 
and a ~ 1. This expression is valid both for systems 
that do and do not undergo a magnetic phase transition. 
In the former case T m is the temperature of the phase 
transition and ST corresponds to the width of the criti- 
cal region around T m (it is such that T m — ST > 0). On 
the other hand, if there is no phase transition (e.g., when 
dealing with the formation of a spin "glass"), one sets 
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T\ = ST = T m = and 9 is the characteristic parameter 
of the ordering. 

The expression (Q) introduced in Ref. f| describes the 
"recovery" effect in a phenomenological way (a micro- 
scopic treatment will be presented elsewhere). As can 
be seen from Eq. (0) the spin-flip scattering amplitude is 
constant for T ^s> 9 (or T m ) and decreases continuously 
near and below T ~ 9 (or T m ). Below T\ the scattering 
amplitude T s is a constant independent of temperature. 
A change of only several percents of the value of T s (T) at 
low temperatures can strongly affect the temperature de- 
pendencies of properties such as H C 2 or J C (T) (see Ref. || 
and below). 

Generally speaking, the theory contains three param- 
eters: To, (3 and 9. To determine their value one has to 
consider separateley the case of conventional and high-T c 
superconductors. In the first case, the only existing ex- 
periments studying the effect of ordering on pair-breaking 
were done on systems that are (anti)ferromagnetically 
ordered. A reasonable value for the characteristic tem- 
perature 9 is thus given by the critical temperature of 
the phase transition which is of the order of IK. Fur- 
thermore, the value of T S (T C ) can be obtained within the 
Abrikosov-Gor'kov theory, once the experimental value 
of the critical temperature T c in the presence of mag- 
netic impurities is given. Its value corresponds to the 
high-temperature limit of the spin-flip scattering ampli- 
tude (when 9 <C T c ). We are thus left with one free 
parameter. 

In the case of high-temperature superconductors the 
analysis of H c i performed in Ref. || delivers the values 
for all parameters. There is thus no free parameter left 
for the high-Tc materials studied in Ref. [|. For exam- 
ple, for overdoped T^E^CuOg studied experimentally 
in Ref. || one has r = 95K, (3 = 1.26, 9 = IK, whereas 
for Bi 2 Sr 2 Cu0 6 studied in Ref. |l| one has T = 105K, 
j3 = 1.38 and 9 = 1.6K. 

Based on Eqs. (jl|)-(|) and (§ we can calculate the am- 
plitude of the Joscphson critical current in the presence 
of magnetic-impurity ordering at low temperatures. 



III. RESULTS AND DISCUSSION 



We consider the weak and strong coupling cases sepa- 
rately. Since, to our knowledge, no experimental study of 
the ordering trend on J C (T) has been undertaken either 
on conventional or high-temperature superconductors for 
systems that do not undergo a magnetic phase transis- 
tion, we chose arbitrary but realistic values of the param- 
eters and discuss qualitatively what can be expected for 
the two types of superconductors. 



A. Weak Coupling 

For a weak-coupling superconductor Z = 1, A(iuj n ) = 
A and one neglects the dependence on u> n i of D nr u.. 
Eqs. (0)-(||) reduce to the Abrikosov-Gor'kov theorytl. 
For the numerical calculations we determine A and T s 
from the experimental values of T c and T c o (the criti- 
cal temperature with and without magnetic impurities). 
A is obtained from Eqs. (0)-(|) for T s = 0, i.e. A = 
In -1 {tt/T c0 ). On the other hand, T S (T C ) can be ex- 
tracted from the Abrikosov-Gor'kcw equation for T c in 
the presence of magnetic impuritiesu: 




where 7 S = r s (T c )/27rT c is the pair-breaking parameter 
and ip is the psi function. 

Fig. 1 shows the temperature dependence of the 
Josephson amplitude for different concentrations of mag- 
netic impurities, that is for different values of Tq. The 
temperature and amplitude of the current are normalized 
to the critical temperature T c o and the amplitude J^ B (0) 
[see Eq. ([)])] in the absence of magnetic impurities, re- 
spectively. As seen in the figure, the presence of magnetic 
impurities lowers the value of T c and J C (T) with respect 
to the Ambegaokar-Baratoff result. As stated in the in- 
troduction this is a consequence of the pair-breaking ef- 
fect. One notes further that J C (T) does not saturate as 
T — > but increase in a temperature range given by 9. 
The upturn is the direct signature of the impurity-spin 
ordering trend, that is, the temperature dependency of 
the spin- flip scattering amplitude. There is a "recovery" 
effect due to impurity-spin ordering opposing the satu- 
ration trend of J c for uncorrelated spins. Note that the 
recovery is not complete [J C (T = 0) < J^ B (0)] since 
T S (T — > 0) 7^ 0. This is even true for a fully ordered 
state. The mechanisms impeding the total recovery of 
the Ambegaokar-Baratoff result (as, e.g., dipole-dipole 
interactions; see also Ref. |l6|) will be described elswhere. 

The ordering trend of local moments leads to a posi- 
tive curvature of J C (T) at low temperatures. Fig. 1 shows 
that, the high-temperature linear part becomes dominant 
as one increases the concentration of magnetic impurities. 
The account of correlations leads thus to a behaviour of 
J C (T) that is qualitatively different from the result ob- 
tained for H C 2(Tp. Indeed, in the absence of spin cor- 
relations, J c and H C 2, considered as functions of temper- 
ature, both have a negative curvature up to T c and are 
linear in the vicinity of T c . Adding magnetic impurities 
extends the high-temperature linear part of J C (T/T C ) to 
lower temperatures to the expens of the low-temperature 
region (that has positive curvature if impurity correla- 
tions are taken into account; see Fig. 1). The effect on 
H C 2(T) is opposite. For high enough magnetic impurity 
concentrations, H C 2(T) displays a positive curvature over 
the whole temperature range when correlations are taken 
into accounta The concentration of magnetic impurities 
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at which a positive curvature of H C 2 was observed exper- 
imentally corresponds to the case To = 8 of Fig. 1 (i.e. to 
the gapless regime). 

Fig. 2 gives another representation of the situation pre- 
sented in Fig. 1. Compare first the solid and dashed lines. 
The dashed curve corresponds to the case Fo/T c o = 0.6, 
6/T c o = 0.1 and (3 = 1.3 (same parameters as the dot- 
ted line of figure 1) whereas the solid line corresponds to 
the case of uncorrelated spins and r s /T c o ~ 0.73. The 
latter value was chosen in such a way that the critical 
temperature T c is the same for both curves. One can 
see that the result obtained by taking into account the 
correlation of impurity spins strongly deviates from the 
curve expected when the spins are uncorrelated. This 
deviation is already significant at temperatures near T c 
for high enough impurity concentrations, but is largest 
near T = 0. Furthermore, the amplitude of the deviation 
is such that it should be experimentally observable. 

Another conclusion can be drawn when comparing the 
dashed and dotted lines of Fig. 2. These curves were ob- 
tained for different values of the parameter (3 [see Eq. (0)] . 
Thus, both curves include the spin-correlations of the 
magnetic impurities, but the strength of the correlations 
is different. The dashed line corresponds to f3 = 1.3 
whereas the dotted line is for (3 = 1.5. Because Tq re- 
mained unchanged, the two curves have slightly different 
T c 's. The increase of (3, which corresponds to an increase 
of the influence of magnetic-impurity correlations, im- 
plies a stronger temperature dependency of J C (T). Fig. 2 
shows that impurity correlations reduce the high and 
intermediaire temperature negative curvature of J C (T). 
Note, finally, that the value at T = is independent of 
(3. 

It would be interesting to study experimentally the 
same high-temperature superconductors as those for 
which H C 2 was measured and described in Ref. [|. In- 
deed, the values of the parameter for the expression of 
r s (T) are given by the analysis of H C 2(T), leaving no 
free parameter (see previous section). Of special inter- 
est would be the case of YBa2(Cui_ a; Zn :E )307_5 since 
there are evidences that the depression of T c and other 
superconducting properties is mainly due to the increase 
in magnetic impurity concentration (see also section on 
strong-coupling superconductors) . 

Note, finally, that our calculations suggest a possible 
application of Josephson junctions based on supercon- 
ductors doped with magnetic impurities (as e.g. high-T c 
materials). Indeed, because the temperature dependence 
of the Josephson current is linear for superconductors 
in the gapless regime (see dash-dotted line of Fig. 1), 
the measure |Of J C (T) for such systems can be used as a 
thermometerlill. 



B. Ternary compounds 

Up to now we have assumed that the magnetic mo- 
ments are correlated but remain disordered. In the case 
of a magnetic phase transition, the result differs from 
above in that the temperature dependency of T s van- 
ishes below T m [see Eq. (0)]. This is shown on Fig. 3 
for parameters corresponding to the experimental results 
on SmRh4B 4 Ll. The critical temperature in the absence 
of magnetic impurities is set to T c0 — 11.4K, the value 
of T c o for LURI14B4 which has the same structure as the 
Sm compound, is also superconducting but has no mag- 
netic moments. The values T c = 1.815K, T m =p0.87K, 
ST = 0.21K were taken from the experiment!! The 
value Tq — 9.75K was determined from the measure of 
J c (Ti = T m — 8T). Finally, the remaining parameter 
(3 = 1.029 was fixed by the value T S (T C ) chosen in such a 
way that T c ~ 1.815K from Eqs. fjj) and (§). The values 
of all parameters have thus been taken from the experi- 
ment and there is no free fitting parameter left. Regard- 
ing the determination of these parameters, one notes first 
that we are in the gapless regime (f = V S (T) / A(T) > 1) 
and thus in the situation of the dash-dotted curve of 
Fig. 1. In addition, the value of (3 that arises because of 
the spin-spin correlations is close to one (value at which 
the spins are uncorrelated) and thus much smaller than 
most of the high-T c superconductors studied in Ref. |] 
(which are of the order (3 ~ 1.3). 

Comparing our result with the experiment on the 
ternary compound, one notes that the behaviour above 
the temperature T± (at which the transition to the an- 
tiferromagnetic state is complete) is well reproduced by 
our calculation. This demonstrates that the increase of 
the current observed near T\ in Ref. ^ can indeed be ex- 
plained by the ordering trend of the magnetic impurities. 
Below Ti, however, the situation is less clear. If we use 
Eq. (g) then we observe a distinct discrepancy between 
experiment and theory. In the experiment, the current 
continues to increase as T is lowered, whereas the cur- 
rent saturates in our calculation (see dashed line below 
Ti in Fig. 3). One reason for this discrepancy lies in the 
fact that in metallic systems containing magnetic impu- 
rities several processes occur that were not accounted for 
in the present analysis. One of them is connected with 
symmetry breaking caused by the impurityQ Another 
reason for the discrepancy is related to the direct dipolc- 
dipole or RKKY type interaction between the impurity 
spins. In the specific case studied in Ref. [?], the occurence 
of a proximity effect due to unoxidized parts of the Lu 
layer of the junction may also be responsible for the lin- 
ear behaviour of J C (T) for T < T x (see Ref. 0). Further 
experimental studies should be carried out on this sys- 
tem to determine which mechanism is the most relevant 
to explain the data. 

To take into account the additional temperature de- 
pendence of the spin-spin correlations induced by the 
above processes, we have considered a temperature de- 
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pendent T s below T\ as well, but with modified values 
of the parameters. The experimental value J c (Ti) gives 
r s (Ti), 9 is kept unchanged and (3 = 1.08 gives the best 
fit to the data. As seen on Fig. 3 (solid line), the temper- 
ature dependence of T s below T\ can account well for the 
experimental result. This is not a trivial statement be- 
cause, as shown in the previous figures, the temperature 
dependency J C (T) cannot be modified at ease by chang- 
ing the values of the parameter. The fact that one has to 
include a temperature dependency of T s in the antiferro- 
magnetic region leads us to the conclusion that not only 
the contact interaction but also the other mechanisms 



(see above and Rcf. 16) are important in this material. 



C. Strong Coupling 

We study the influence of strong coupling onto the re- 
sults of the last section. This case is numerically more 
difficult to handle, due to the fact that the renormal- 
ization function has to be considered and that all three 
unknown functions A„, Z n and u n depend on n. As be- 
fore, we determine T s and A from the equations for T c and 
T c0 . For T c one linearizes Eq. (@). As for T c0 (T s ^ 0) 
we use the expression derived by one of the authorala: 
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where fi is the characteristic phonon frequency. Fig. 4 
shows a generic case for the strong coupling limit. There 
is no qualitative difference between the weak and strong 
coupling results. The parameters in Fig. 4 were chosen so 
as to correspond to the case^if YBa^Cui^Zn^aOy-s 
with x ~ 0.04 (T c ~ 55K)B Q . This result leads us 
to conclude that the overall qualitative behaviour ob- 
tained for weak-coupling superconductors is not modi- 
fied by increasing the strength of the coupling. In par- 
ticular, one expects a nearly linear behaviour of J C (T) 
for high-T c materials where T c has been strongly reduced 
by the presence of magnetic impurities and are in the 
gapless regimCr-Examples of such systems, are Zn-doped 
YBa 2 Cu 3 7 -Jia a and Lai. 85 Sro.i 5 Cu0 4 lia 6 > c , as well as 
other under- and overdoped materials studied in Ref. |] 
(Tl 2 Ba 2 Cu0 6 , Bi 2 Sr 2 Cu0 6 , or Sm 1 . 85 Ce . 15 CuO4- !/ ). It 
would be interesting to perform measurements on these 
(or other) systems to verify our predictions. As stated 
earlier, the above mentioned materials are of special in- 
terest, since all parameters necessary to calculate r s (T) 
and J C (T) [Eqs. (0) and ([|l respectively] were deter- 
mined from the study of H^Q- 



IV. CONCLUSION 

In summary, we have calculated the amplitude of the 
Josephson current J C (T) in the presence of magnetic im- 
purities taking into account their ordering trend at low 



temperatures. We have shown that the "recovery" effect 
introduced in Ref. || is also occuring for the Josephson 
current J C (T). The effect manifests itself in an upturn 
of J C (T) at low temperatures (instead of the saturation 
occuring for uncorrelated magnetic moments). The effect 
is sizeable and qualitatively similar in weak and strong- 
coupling superconductors. The theory was used to ex- 
plain the experimental results obtained for SmRli4B4 and 
leads to good agreement for temperatures above the mag- 
netic phase transition. We have also presented results for 
high-T c materials. 

It would be interesting to measure the effect in both 
conventional and high-temperature superconductors so 
as to verify our predictions. The effect might be larger 
in the latter (for overdoped or Zn-doped systems) as ii*- 
dicated by the unusal temperature dependence of if c2 Q, 
since most of these materials are in the gapless regime. 
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Figure captions: 



Fig. 1: Temperature dependence of the Josephson current amplitude J C {T) for different magnetic impurity con- 
centrations. The temperature and the current amplitude are normalized to the values T c q and J^ B (T = 0) in the 
absence of magnetic impurities, respectively. (3 = 1.3, 9/T c q = 0.1 and Tq/T c q — (Ambegaokar-Baratoff; solid line), 
r /T c0 = 0.4 (dashed), r /T c0 = 0.6 (dotted) r /T c0 = 0.8 (dotted-dashed). 

Fig. 2: Temperature dependence of the Josephson current for different values of (3. Solid line: (3=1, dashed line: 
= 1.3, dotted line (3 = 1.5. r /T c0 = 0.6, 9/T c0 = 0.1. 

Fig. 3: Josephson current for SmRl^Bzi. Solid line: theoretical result with T S (T), dashed line: J C {T) for 
r s (T < Ti) = To, dots: experimental points from Ref. ^ The parameters are given in the text. 

Fig. 4: Influence of impurity-spin correlations on the Josephson current amplitude for strong-coupling superconduc- 
tors, solid line: uncorrelated spins, dashed line: correlated impurity spins. The parameters correspond to the case of 
YBa2(Cui_ 2; Zn :I ;)307_5 with x ~ 0.04 (see text). 
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